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£S) _ Abstract. The first 3D calculation of shock wave propagation in a homogeneous 

■ QGP has been performed within the new formulation of relativistic dissipative 

, hydrodynamics which preserves the causality. We found that the relaxation time plays 

an important role and also affects the angle of Mach cone. 
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1. Introduction 

One of the most important questions to be clarified in the hydrodynamical approach 
to the relativistic heavy ion physics is the effect of dissipative processes. The second 
i-£h ' order theory, first proposed by Muller and developed by Israel and Stewart, has been 
considered standard approach for this problem pQ . But it is quite complex and involves 
^1 many unknown parameters from the point of QCD dynamics so that its complete 

application to practical problems such as relativistic heavy ion reactions has not been 
done yet [2]. 

In this work, we propose an alternative approach to this question [3]. We show 
that the physical origin of the second order theories can easily be understood in terms 
of memory effects. The irreversible current modified by the memory effects becomes 
consistently with causality and sum rules pE]. Based on this idea, we introduced the 
memory effect to the relativistic dissipative hydrodynamic of Landau [5], where we 
introduce only one extra parameter, the relaxation time T re i ax , in addition to the usual 
viscosity coefficients of the Navier-Stokes equation. The resulting equation becomes 
hyperbolic [3]. 

The effect of viscosity is also important when we discuss the possible generation and 
propagation of shock waves in the QCD medium created in the process of relativistic 
heavy ion collisions. As discussed extensively in this conference [6], it has been suggested 
that a high energy jet propagating in the QGP may generate a Mach cone and 
observables associated with such phenomena may bring important information of the 
genuine hydrodynamical properties of the matter [7]. The dynamical simulation of 
shock wave generation is very difficult even for the non-relativistic regime. A full 3D 
simulation of shock wave dynamics has never been done for the heavy ion collisions. In 
this work, we apply our formulation to the calculation of full 3D relativistic (causal) 
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shock wave problem. The implementation of our method to the existing ideal hydro- 
codes is straightforward, particularly to those based on the local Lagrangian coordinate 
system such as SPheRIO [8]. We organize the present work as follows. In the next 
section, we briefly introduce our formalism and discuss its application to the generation 
of shock waves. In Section 3, we present some results of 3D calculation of shock wave 
propagation within the causal formulation of the dissipative hydrodynamics. In Section 
4, we discuss the result and perspectives. 



2. Causal Approach to the Relativistic Dissipative Hydrodynamics 

The fundamental problem of the first order theory like the Navier-Stokes theory comes 
from the fact that the diffusion equation is parabolic. The physical origin of this problem 
can be followed up to the fact that the irreversible current j is assumed to be proportional 
to a thermodynamic force F as 

J = -LF, (1) 

where the Onsager coefficient L is, in general, a function of thermodynamic quantities. 
Usually F is related to the inhomogeneity in the density. When the microscopic 
rearrangement time scale is not negligible compared to the time scale of the change 
in the irreversible current, then the above should be replaced by the equation of motion 
for the current, 

! = --iFM--J. (2) 

9t T R T R 

where r R is the relaxation time. For very small t r , we recover Eq.([T]). Thus Eq.([T]) 
can be understood as the large viscous limit of the damped motion, where the velocity 
(current) is proportional to the force (Aristotelian vision). It can be shown that the 
above modification is enough to convert the parabolic nature of a diffusion equation to 
hyperbolic one [3]. 

For the relativistic hydrodynamics, we have to consider several different kind 
of thermodynamical forces related to the velocity and density inhomogeneity. They 
are F = d a u a , F^ = d^ot and F^ u = d^u V) where is the four- velocity of the 
fluid, a = fi/T with the chemical potential. These inhomgeneities generate the 
corresponding irreversible currents and the analogous equations to Eq.Q for them 
should be 

r/TT dTT^ du^ 

T R — = -Cd a u a -U, t r — = r)dW-^ v , r R — = -«^a-^,(3) 
dr dr dr 

where djdr = u^d^, and (, t] and k are bulk viscosity, shear viscosity and thermal 
conductivity coefficients, respectively. The energy-momentum tensor is expressed with 
these currents as T^ v = eu^u v - (p + IT) P^ + P^ al3 iT a p, where P^ a = g^ u - u^u v and 
p^vp ig the double symmetric traceless projection, P^ afi = \ [P^ a P v < 3 + prfp™} _ 
P Mi/ P Q/3 /3, whereas the conserved baryon number current is given by = nu^ + P^ a v a . 
The hydro dynamic equations are 

<9 M T^ = 0, <yV = 0. (4) 
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3. Shock wave 

Although Eqs.( l3lHl) together with the equation of state give the complete description of 
the hydrodynamical motion of the system, in practice, some additional care to be taken, 
especially for the simulation of the shock wave dynamics. Whenever there exists a shock 
wave, always occurs an entropy production through the shock front. In an idealized 
hydrodynamical approach, the shock front is a discontinuity in thermodynamical 
quantities in a hydrodynamic solution. Mathematically speaking, it should be treated 
as the boundary condition to connect two distinct hydrodynamic solutions. Physically, 
it is not a real discontinuity, but a quick change of the density in the region where the 
local equilibrium is not satisfied. Thus it has a finite thickness at least a few times of the 
mean-free path (typical microscopic scale of distance) for a stationary shock. Under a 
dynamical condition such as relativistic heavy ion collisions, the compression shock may 
have much more larger thickness due to the many complicated local transient properties. 

To reproduce true shock wave phenomena, the full degrees of freedom of the 
hydrodynamics, together with a proper boundary condition correctly connecting to 
regions through the non-equilibrated domain of the shock, are required. The usual 
numerical approach of hydrodynamics excludes such a possibility from the beginning. 
Since there exist no short wavelength excitation modes due to the finite discretization, 
the energy and momentum conservation required by the hydrodynamics result in a 
very rapidly oscillations of the variables near the shock region. These quick oscillations 
are somewhat compensating the thermal energy associated to the entropy production 
throughout the shock front. In order to avoid these oscillations, von Neumann and 
Richtmeyer[9] introduced the method of pseudo- viscosity. The idea is to put the bulk 
viscosity in the shock region to replace the entropy production through the shock front. 
The form of the pseudo-viscosity is chosen so that asymptotic values of thermodynamic 
quantities connects smoothly from one side to the other through the shock, satisfying 
the Huguniot-Rankine boundary condition. For the relativistic extension, we have to 
transform this formalism in the causal form. Thus, we introduce the causal pseudo- 
viscosity by using the bulk viscosity coefficient ( in Eq.Q as ( = — Co s [1 ~ otO] for 
9 < and otherwise where 9 = and (q and a are constants proportional to the 
space resolution scale of the numerical solution. From the condition of causality, the 
relaxation time associated to this pseudo- viscosity should satisfy > 2 ^^ .This means, 
the larger the viscosity, the larger the relaxation time. In Fig.l, we show the example of 
creation and propagation of shock wave induced by a small QGP drop injected into the 
homogeneous QGP fluid at a relativistic energy (7 = 10 ). Here, we take a massless ideal 
gas equation of state, and the incident drop has half the temperature of the background 
in the local rest frame. The relaxation time necessary to maintain the causality in this 
calculation is found to be 5/m/c. This seems to be large, but the effect of shock wave 
requires a large viscosity so that a large relaxation time is necessary. We can see the 
generation of Mach cone and also the turbulent structure after the shock cone. We 
found that the cone angle depends on the viscosity and also the related relaxation time 
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tr. For a finite tr the Mach cone opens more than the usual estimate for ideal fluid 
case. 

4. Discussion and Perspectives 

We presented the first full 3D relativistic viscous hydrodynamic calculation which 
preserve the causality correctly. In the case of shock wave simulation, it is fundamental 
to include the bulk viscosity to take account for the entropy generation through the 
shock front. To satisfy the causal propagation of the signal, the presence of viscosity 
requires a finite relaxation time. It is found that the Mach cone angle depends on the 
viscosity and the relaxation time. It is also observed that a very turbulent structure 
after the shock. More quantitative and systematic study of the shock wave generation 
and its propagation in an expanding QGP is in progress. 

This work has supported by CNPq, FINEP, FAPERJ and CAPES. T.K. appreciates 
the discussion and encouragements of H. Stoecker, E. Fraga and A. Muronga. 
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